Decay of the -0(3770) to Light Hadrons 
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Abstract 

We describe a search for -0(3770) decay to two-body non-DD final states in e + e" data produced 
by the CESR collider and analyzed with the CLEO-c detector. Vector-pseudoscalar production 
p°ir°, p + 7r~, ujit , (f)ir°, pn, un, <pn, pn', ton' , 4>n' , K*°K°, and K* + K~ is studied along with that 
of 6i7r (bi7T° and bfn~) and 7r + 7r~7r°. The largest amount of disagreement between the expected 
rate for e + e~ — ► 7* — > X and that for e + e~ — > X at yfs = 3.773 GeV is found for X = 4>n, at an 
excess cross section of (2.4 ± 0.6) pb [r^ 3770 ^ = (7.4 ± 1.6) keV], and a suggestive suppression is 
seen for 7r + 7r~7r° and pir. We conclude with form factor determinations for lott , pn, and pn'. 
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The ?/>(3770) charmonium state decays most copiously into the OZI-allowed DD pair 
owing to the closeness of the mass threshold. Hadronic or radiative transitions to lower- 
lying charmonium states, decay to lepton pairs, or decay to light hadrons are all available, 
but their branching fractions are highly suppressed. The ■0(25')-'0(377O) mixing scenario 
proposed in [l| gives rise to an enhancement of certain partial decay widths, but the resulting 
branching fractions are still small due to the large width of the -0(3770). Nevertheless, some 
of the branching fractions are within experimental reach, and experimental information on 
0(3770) non-DD decays has recently begun to emerge [2j. 

This Letter describes the search for -^(3770) decay to vector pseudoscalar (VP) final 
states (p°7r°, p + n~, un , prj, ur], <prj, pi]', urf, <fir)', K*°K°, K* + K~) in CLEO-c data taken 
at the ■0(3770) resonance. We also seek 7r + 7r _ 7r° as it is a mode that exhibits curious struc- 
ture in ip{2S) decay and bin (in both the charged and the neutral isospin submodes) as 
the most commonly produced two-body hadronic final state in ip(2S) decays. We use data 
samples taken at two energies, y/s = 3.671 and y/s = 3.773 GeV. We establish event yields 
in both by counting events satisfying the selection criteria detailed below, and subtracting 
misreconstructed and therefore erroneously selected events. We measure the visible cross 
section at both center-of-mass energies for all modes. The sideband-subtracted event counts 
at y/s = 3.773 GeV are compared with the expected rate from e + e~ — > 7* — > VP (contin- 
uum), in order to discern a statistically significant discrepancy between the two. Assuming 
the continuum cross section a(e + e~ — > 7* — > VP) is given by 

„( e \ - 4to cm \F(s)\ 2 4p( S ) m 
a \ b ) ~ 3 s 3/2 ' \ L ) 

its measurement gives access to the form factor T . The momentum of either hadron is 
denoted by gyp- F° r an channels, the event yield at y/s = 3.671 GeV is solely due to the 
above process. Also, for certain channels the event yield at y/s = 3.773 GeV will be entirely 
attributable to continuum production, namely those that cannot be produced through cc — > 
ggg because of isospin suppression, such as um°, prj, and prj'. Their remaining open avenue 
for 0(3770) decay is cc — > 7*, which is severely suppressed. 

We use e + e~ collision data at y/s = 3.773 GeV (£ = 281 pb -1 ) and y/s = 3.671 GeV 
(C = 21 pb _1 ). The data analyzed here were collected with the CLEO detector operating 
at the Cornell Electron Storage Ring (CESR) The CLEO detector features a solid 
angle coverage of 93% for charged and neutral particles. The charged particle tracking 
system operates in a 1.0 T magnetic field along the beam axis and achieves a momentum 
resolution of ~ 0.6% at momenta of lGeV/c. The Csl crystal calorimeter attains photon 
energy resolutions of 2.2% for E 1 = 1 GeV and 5% at 100 MeV. Two particle identification 
systems, one based on energy loss (dE/dx) in the drift chamber and the other a ring- imaging 
Cherenkov (RICH) detector, together are used to separate kaons from pions. The combined 
dE/dx-RlCH particle identification procedure has a pion or kaon efficiency >90% and a 
probability of pions faking kaons (or vice versa) <5%. 

We identify intermediate states through the following decays: n° — > 77, 77 — >• 77 or 
7T + 7r~7r°, T]' — > n + n^ri (rj — > 77 only), K° Ks — > cu — ► p — > , 

<f) — > K + K~ or 7r + 7r~7T°, K* — > Kw, and h\ — > ujtt. Event selection proceeds exactly as for 
CLEO's ip(2S) -> h x h 2 analysis 0, the main features of which were requirements on the 
total energy E y { s and momentum as well as on the invariant masses of intermediate particles, 
in combination with particle identification criteria. For some modes which are particularly 
susceptible to background from radiative Bhabha or p + p~ production, we tighten the se- 
lection in the present analysis by imposing the following additional requirements: For p + 7i~ 
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and 7r + 7r~7r°, e + e~ — > events with a fake n° candidate are suppressed by a decay 

angle requirement of | cosa| < 0.8. For pq, backgrounds from e + e~(7) final states with a 
fake 7] candidate are reduced by allowing neither 7r =1= to satisfy electron identification criteria. 

We present distributions of scaled total energy E vis / y/s and reconstructed invariant 
masses for selected modes in Figures All cuts other than that imposed on the quantity 
displayed have been applied. 

The efficiency e for each final state is obtained from signal Monte Carlo with the 
EvtGen event generator 0], including final state radiation [7], and a GEANT-based detec- 
tor simulation We generated the VP modes with angular distribution (1 + cos 2 #) 0, 
bin flat in cos#, and 7r + 7r~7r° as in u decay. We assume B(b\ — > co>7r)=100%. 

Systematic uncertainties on the cross section measurements arise from various sources, 
some common to all channels, some channel specific: The systematic errors on branching 
fraction ratios share common contributions from the uncertainty in luminosity (1%), trigger 
efficiency (1%), and electron veto (0.5%). Other sources vary by channel, including cross- 
feed adjustments (50% of each subtraction), MC statistics, accuracy of MC-generated polar 
angle and mass distributions (5% for bin, 14% for 7r + 7r~7r ), and detector performance 
modeling quality: charged particle tracking (1%/track), n°/rj and K s finding (2% / (tt / rj) , 
5%/Ks), ir/K identification (3%/identified ir/K), and resolutions of mass (2%) and total 
energy (1%). 

Systematic uncertainties dominate in the cross section measurements for most channels 
at y/s = 3.773 GeV data and are comparable to the statistical errors for some modes at 
y/s = 3.671 GeV. 

The signal yields at both center-of-mass energies are listed in Table EJ separated into 
signal mass windows and sideband counts. Also listed are the efficiencies and cross sections. 
The statistical errors arise from 68% CL intervals. All cross sections include an upward 
correction of (20 ± 7)% to account for initial and final state radiation effects |3j. In case of 
the isospin- violating modes, we also correct for electromagnetic interference between the tails 
of the J/ip, i/j(2S), and "0(3770) resonances with continuum production by a 4.9% upward 
[1.2% downward] adjustment to the cross-sections at y/s = 3.671 GeV [y/s = 3.773 GeV] |l0(. 
The results in Table H] for y/s = 3.671 GeV supersede those in [3J]. 

We now focus on the discrepancy between the y/s = 3.773 GeV yield and expected con- 
tinuum contribution in order to determine whether there is significant production from 
■0(3770) decays. To arrive at an estimate for the continuum background at y/s = 3.773 GeV, 
two routes are pursued: Method I. We scale the measured yield (after sideband subtraction) 
at y/s = 3.671 GeV by the luminosity ratio, the ratio of efficiencies (0.88 — 1.00), and an 
assumed dependence of 1/s 3 of the continuum cross section, corresponding to a form fac- 
tor dependence of 1/s. This method uses data as much as possible, but suffers from the 
limiting event yield in the lower energy data sample. Using a different power l/s n results 
in a relative change of 5.4% in the scale factor per unit of n. Method II: We use a SU(3)- 
based scaling prediction, whereby the the cross sections <r(e + e~ — > VP) are linked ^l| as 
uin \ pTj '. K*°K° : pn : pq 1 : <jyq : K* + R- : (fyrjf : ujt] : urf : <M = 1 : 2/3 : 4/9 : 1/3 : 
1/3:4/27: 1/9:2/27:2/27: 1/27 : 0. By combining our data of the two isospin violating 
modes with highest statistics, un° and pr\ (scaled up by a factor of 3/2), we determine a 
unit of cross section as a su ^ = (15.1 ±0.5) pb at y/s = 3.773 GeV. This results in a precise 
prediction for each channel, albeit a model-dependent one. We note satisfactory agreement 
between the yields expected on this basis and those observed in the data at y/s = 3.671 GeV 
for all channels except of K*K (also see [3])- No such prediction is made for and bin. 
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For each channel, both continuum predictions are compared with the yield at ^/s = 
3.773 GeV by a method similar to that proposed in l3| • It is the same procedure that was ap- 
plied to i/j(2S) decays in [3]: The probability that the the continuum production from either 
method together with the misreconstruction background as estimated from the sidebands 
fluctuate to an event count equal to or beyond the observed signal yield at t/s = 3.773 GeV 
is calculated with simulated trials governed by Poisson statistics. These, expressed in units 
of standard deviations, are included as S l and S 11 in Table H] We find statistical agreement 
between the yields, with a few exceptions. The mode <firj is found to be enhanced over the 
prediction from either method: The weighted mean excess over continuum production is 
(61.6 ±11.7) events. This corresponds to a cross section of cr^ 3770 -* = (2.4 ±0.5 ±0.3) pb, or, 
using <7 obs (?/>(3770) — > DD) = (6.39 ±0.20) nb ljj] and removing the radiative correction fac- 
tor in a^ 3770 \ a branching fraction #(-0(3770) — > (j)r]) = (3.1 ± 0.6 ± 0.3 ± 0.1) x 10 -4 , where 
the first error is statistical, the second systematic arising from this measurement, and the 
third that induced by a^% 770) . A partial width of rJJ 3770) = (7.4±1.6) keV follows Q. Sug- 
gestive suppressions are observed for 7r + 7r~7r°and p°ir°. The observed K*°K° cross section at 
y/s = 3.773 GeV is consistent with being saturated by the expected continuum production 
as extrapolated from yj~s = 3.671 GeV (Method I). As the observed K*°K° cross sections 
at both energies jj3] far exceed the SU(3) predictions (Method II), and by similar amounts, 
there is no indication of a substantial 0(3770) —*K*°K° contribution. Rather, the excesses 
originate in the continuum process, e + e~ —>K*°K°. The same comments apply to K* + K~, 
but with respect to an observed deficit obtained with Method II. 

Additional information on 7r + 7r~7r° is shown in Figure 0] The dipion invariant masses in 
y/s = 3.773 GeV data shows features similar to that of y/s = 3.671 GeV (i.e., population 
of the p mass bands together with an accumulation at higher masses); the yield reduction 
appears uniform in the dipion invariant mass distribution. 

We compute upper limits on the event yields originating from 0(3770) decays for all 
modes, where we treat those with a deficit as zero counts, neglecting interference effects, 
and arrive at upper limits on the observable cross section excess over continuum as included 
in Table |U 

The measured cross sections for c<j7r , pi], and pr( are converted into form factor measure- 
ments, which are listed in Table |H] Our results are in agreement with, but more precise 
than, those recently reported by BES (ljj . 

In summary, we have sought twelve vector pseudoscalar final states in data at \J~s = 
3.773 GeV. Combined with data collected at y/s = 3.671 GeV, we establish cross section 
measurements for these channels at both energies. We find evidence for the decay -0(3770) — > 
4>r], and see hints that -0(3770) decays to pn and 7r + 7r~7r could be causing a deficit to appear 
in their yields through negative interference with continuum production 0,0]. Otherwise, 
we note broad agreement with the continuum predictions. Form factor measurements for 
ujtt , pr], and prj' have been presented. All our measurements are either firsts of their kind 
or constitute an improvement over previous measurements. 
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TABLE I: The number of events N in the mass signal windows ("sw") and sidebands ("sb") 
in data taken at \fs = 3.671 GeV and y/s = 3.773 GeV data; the efficiency e in percent; the 
level of consistency or significance, expressed in units of standard deviations, between continuum 
background and observed yield, for the two methods of determining the continuum background 
described in the text, S l and S 11 ; the cross sections at \fs = 3.671 GeV and yfs = 3.773 GeV; the 
cross section ^(3770) — > h±h2, computed as the excess over the continuum prediction as established 
using Method I or Method II (see text). 
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TABLE II: Form factors with statistical and systematic errors. 
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FIG. 1: Scaled visible energy E v i B /\/s for selected final states. Circles: data at y/s = 3.773 GeV, 
shaded histogram: data at y/s = 3.671 GeV scaled by luminosity, dashed histogram: signal MC, 
arbitrary normalization. Arrows indicate selection intervals. 
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FIG. 2: Reconstructed invariant mass distributions for selected final states. Symbols as in Fig. 1. 
The figures on the left (right) side pertain to the first (second) of the two final state particles. 
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FIG. 3: Mass distributions for selected final states, continued. Symbols as in Fig. 1. 
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FIG. 4: Dipion invariant mass distributions for the 7r + 7r~7r° final state in data (a) at yfs = 
3.773 GeV, (b) at ->/s = 3.671 GeV. (c) The invariant mass of all pion pairs per event and (d) the 
reconstructed 7T° mass, in data at yfs = 3.773 GeV (circles), data at -y/s = 3.671 GeV, scaled by 
luminosity (shaded histogram), and phase space MC (dashed line). 
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